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ABSTRACT 



Context. The standard cooling models of neutron stars predict temperatures of T < 10* K for ages f > 10' yr However, the likely 
thermal emission detected from the millisecond pulsar J0437-4715, of spin-down age ?, ~ 7x 10' yr, implies a temperature T ~ 10' K. 
Thus, a heating mechanism needs to be added to the cooling models in order to obtain agreement between theory and observation. 
Aims. Several internal heating mechanisms could be operating in neutron stars, such as magnetic field decay, dark matter accretion, 
crust cracking, superfluid vortex creep, and non-equilibrium reactions ("rotochemical heating"). We study these mechanisms to es- 
tablish which could be the dominant source of thermal emission from old pulsars. 

Methods. We show by simple estimates that magnetic field decay, dark matter accretion, and crust cracking are unlikely to have a 
significant heating effect on old neutron stars. The thermal evolution for the other mechanisms is computed with the code of Fernandez 
and Reisenegger. Given the dependence of the heating mechanisms on the spin-down parameters, we study the thermal evolution for 
two types of pulsars: young, slowly rotating "classical" pulsars and old, fast rotating millisecond pulsars. 

Results. We find that magnetic field decay, dark matter accretion, and crust cracking do not produce any detectable heating of old 
pulsars. Rotochemical heating and vortex creep can be important both for classical pulsars and millisecond pulsars. More restrictive 
upper limits on the surface temperatures of classical pulsars could rule out vortex creep as the main source of thermal emission. 
Rotochemical heating in classical pulsars is driven by the chemical imbalance built up during their early spin-down, and is therefore 
strongly sensitive to their initial rotation period. 
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1. Introduction 



- dense matter — stars: rotation — pulsars: general — pulsars: individual: PSR J0437-4715 — pulsars: 



Neutron stars (NSs) are compact objects composed of a liquid 
core enveloped by a solid crust. The core is expected to contain 
superfluid neutrons and superconducting protons, while the crust 
contains heavy atomic nuclei arranged in a crystal lattice, coex- 
isting with superfluid neutrons in its inner part. The high density 
of the NS core, up to ~ (3 - 9)po, where po is the saturation 
nuclear matter density, cannot be reproduced in terrestrial labo- 
ratories. This turns the NSs into natural laboratories. The study 
[ of their thermal evolution, confronting theory and observation, 
■ provides a useful test for the understanding of the properties of 
' matter at supernuclear density. 

For all standard cooling models dYakovlev & PethicHl2004l) . 
^ neutron stars cool down to surface temperatures Tj < lO'* K 
within less than 10^ yr. Nevertheless, the observation of ul- 
traviolet thermal emission from millisecond pulsar J0437-4715 
(iKarg altsev et al. 2004), whose spin-down age, connected to the 
latest distance of 157 pc dPeller et alj|2008i) . is r^^ ~ 7 x 10*^ yr 
dvan Straten et alj|2001h . shows a surface temperature of about 
~ 10^ K for this pulsar. Hence, a heating mechanism needs be 
to added to the standard cooling models to obtain agreement be- 
tween theory and observation. 

There are several heating mechanisms that can be 
present during the late stages of the thermal evolution. 
These include the frictional m otion of superfluid neu- 
tron vortices dAlpar et alj Il984t [Shibazaki & Lamb 1989; 
iLarson & Linkl[l999h . rotochemical heating (iReiseneggeriil995i. 



1997; Fe rnandez & Reis enegger 2005; Petrovich & Reisenegge 
2010), magnetic field decay (G oldreich & Reisenegger 199. 
Thompson & Duncan 1996; Pon s et alJ 2007 ). and crust crack- 
ing (Ba vm&Pinesi il971i; iChenget alj Il992h . Other mecha- 
nisms, based on more speculative hypo theses, such as a time 
variation of the gravitation al constan t (Jofreetal. 2006), the 
decay of exotic particles (iHannestad et a l. 2002), or the ac- 
cretion of dark matter parti cles dde Lavall az & Fair bairnl 120101 ; 
iKouvaris & Ti nvakoV '2010*). could in principle als o heat old 
neutron stars. In Schaab et al. ( 1993) and lLarson & L ink ( 1999), 
several of the internal heating mechanisms cited are studied 
and confronted with observational data of neutron stars (sur- 
face temperature upper limits). However, the oversimplified de- 
scription of some heating mechanisms (rotochemical heating in 
Schaab et al. 1999) and the until then non-detection of thermal 
emission from neutron stars older than ~ 10* yr, made it impos- 
sible to obtain reliable conclusions. 

The goal of this work is to provide a comparative analy- 
sis of the thermal evolution including different heating mecha- 
nisms. In order to do this, we discard some of them as not strong 
enough (magnetic field decay, dark matter accretion, and crust 
cracking), and we present a more detailed study of the most 
promising ones: vortex creep and rotochemical heating. We con- 
front these mechanisms with the thermal emission detected in 
the millisecond pulsar J0437-4715 and the best available upper 
limits on the temperature of six other old pulsars. Owing to the 
dependence of heating mechanisms on spin-down parameters, 
which leads to different temperatures for different pulsars, we 
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separately study the thermal evolution for two types of pulsars: 
young, strongly magnetized, and slowly rotating "classical" pul- 
sars, and old, weakly magnetic, and fast rotating millisecond pul- 
sars (MSPs). 

The paper is organized as follows. In Sect. |2] we describe 
the different heating processes and assess their importance. In 
Sect. [3] we show the effects of vortex creep and rotochemical 
heating in classical and millisecond pulsars. In Sect. |4]the pre- 
dictions are confronted with observed data. A summary of our 
main conclusions is given in Sect. |5] 

2. Mechanisms 

The cooling of a NS is caused by neutrino emission from the 
interior and by thermal photon emission from the surface of the 
star. The neutrino emission is determined by the state of matter 
in the core, which depends on the stellar mass and the properties 
of the matter at nuclear density. The evolution of the internal 
temperature of an NS is given by the thermal balance equation 



2.2. Dark matter accretion 



(Lh - Ly - Ly), 



(1) 



where C is the total heat capacity of the star, Ly is the photon 
luminosity, Ly is the neutrino luminosity, and Lh is the power 
generated by internal heating mechanisms, of which we consider 
the following. 



2.1. Magnetic field decay 



iGoIdreich & Reisenegger! (Il992h studied the processes that 
promote the dissipat i on of magnetic energy in NSs. Later, 
[Thompson & DuncanI (Il996h studied the resulting emission (X- 
rays, neutrinos, Alfven waves) from very strongly magnetized 
NSs ("magnetars"), and'Po ns et al.l (l2007h showed observational 
evidence that suggests heating by this mechanism in different 
classes of relatively young, strongly magnetic NSs. Because the 
physical processes involved in the decay of NS magnetic fields 
are still uncertain, we only make an order-of-magnitude estimate 
of the magnetic field B required to produce a detectable surface 
temperature in old NSs. 

At ages f > 10^ yr, cooling is dominated by photon emission. 
The luminosity caused by the decay of the magnetic field in a NS 
with radius R and magnetic energy Eg ~ {4-nR^/3){B-}/3:7T in a 
time scale t, is 



L = AnR^^rTt ~ — 



AnR^ (B^) 1 



(2) 



In this way, the magnetic field required to account for an NS with 
a surface temperature T^ ~ 10^ K (i.e. the detected temperature 
in J0437-4715) and radius /? = 10 km is 



Brms = V<B2> 



llA-no-T^t 



1/2 



= lO'^VfyG, 



(3) 



where tj is the age of the NS in units of 10^ yr An old classical 
pulsar of ~ 10^ yr and an MSP of ~ 10^ yr require, respectively, 
magnetic fields of ~ 10'^ G and ~ lO^'* G in order to obtain a 
detectable temperaturqj. Thus, the magnetic fields inferred from 
the spin-down in classical pulsars, ~ 10" G, and MSPs, ~ 10^ G, 
are much too low to produce detectable heating. 



' These values might be reduced if a substantial fraction of the pro- 
tons is the star are superconducting JEasson & Pethic]dll977h . 



Recently, 



Kouvaris & Tinvakovl 



(l20Toh 



and 



Ide Lavaliaz & FairbairnI (1201 Ol) studied the effects of dark 



matter (DM) accretion and annihilation in compact objects. One 
of these is the potential increase of the surface temperature of 
old NSs in sufficiently dense environments (e.g., the Galactic 
center and globular clusters). Because we are interested in the 
pulsars of the solar neighborhood, where the DM density is 
relatively low, we only make an order-of-magnitude estimate 
of the maximum possible thermal emission caused by this 
mechanism. The maximum possible accretion rate onto a NS is 
given by 



M » PDMVcoTTb^ 



(4) 



where poM is the DM density, Voo is the relative velocity between 
DM particles and the NS, and bcc is the maximum impact pa- 
rameter for a DM particle to hit NS. This is valid as long as 
the DM-baryon (or DM-lepton) cross section {cr} > nR^/Nhar ~ 
10 "^^ cm-, so any DM particle entering the NS stays inside. If 
(cr) is smaller than this, M is reduced correspondingly. From the 
Milky Way (MW) rotation curve, the DM density in the solar 
neighborhood Pom ~ 10"^ Mq pc"-^. From Newtonian energy 
and angular momentum conservation, the maximum impact pa- 
rameter is bcc ~ Rvesc/voo, in the limit Vt,sc s> Vco. Considering 
the escape velocity v^sc ~ 2c/3 and the typical velocities for NSs 
moving through the MW, Voo <: 2 x 10^ cm s"', we can obtain an 
upper limit on the DM accretion rate, M ^ 9 x 10"^^ Mq yr '. 
Hence, the maximum luminosity is MGM/R ^ 10^^ erg s"' for 
stable DM , and Mc~ » 5x10^^ erg s"' for decaying DM. 
Considering that the cooling is dominated by photon emission, 
these correspond to the surface temperatures ~ 2 x lO"* K and 
~ 3 X 10^ K for stable and unstable matter, respectively. In order 
to get a detectable T^ ~ 10^ K, a lO'' times higher DM den- 
sity is required (as perhaps found very near the Galactic center). 
Thus DM accretion does not have a significant effect in the solar 
neighborhood. 

2.3. Crust cracking 

This mechanism considers a rotating NS whose crust solidifies 
with an ellipsoidal form. The spin-down causes a gradual change 
to a more spherical shape. In this process, the stress in the crust 
of the NS is increased. However, the rigidity of the crust causes 
it to remain more oblate than it would be if it had no resistance 
to shear. When the crust reaches a critical deformation, it breaks, 
part of the accumulated strain energy is released, and the excess 
oblateness, due to the crust rigidity, is reduced. The mean stress 
Q- in the c rust caused by the spin-down is given by cr = yu(eo - e) 
dBavm & Pines 1971), with e - (I - Io)/Io, where jj. is the mean 
shear modulus of the crust, / is the moment of inertia of the 
star, Iq is its non-rotating value, and eo is a value of e at which 
the crust is stress-free. As the stellar rotation slows down, e is 
correspondingly decreased, increasing the accumulated stress. 
Eventually, the crust cracks, eo is suddenly decreased, causing a 
discrete change A(en - e). The amoun t of strain energy released 
by this crack is JBavm & Pineslll97ll) 



A£,„ = -2B(eo - e) ■ A(eo - e), 

and the time between successive quakes is 



At ^ 2A- 



Q.Q. 



de 



(5) 



(6) 
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where Q. is the angular velocity of the star, B ~ jj^Vd 42 (with Vc 
the crust volume), and A quantifies the increase in gravitational 
energy due to the deformation o f the shape (for more detail, see 
ICutler et alJ2003HZdunik et al.ll20 08). If eo - e =* 6^, the critical 
breaking strain angle, and the time between quakes is small com- 
pared to the timescale of thermal evolution, the time-averaged 
energy dissipation rate by this process is 



Lcc = biOcniCii, 



(7) 



', where P_2{) is the period derivative mea- 
^^ and Psms is the period in units of 5 mil- 



where b = B/A. Thi s expression, based on the formalism o f 
iBavm & Pinesl (Il97lh . differs from tha t of ICheng et all (Il992h . 
but agrees with the correction made by 'Schaa b et alj (11999') for 
the particular case of a neutron star modeled as a Maclaurin 

sp heroid. 

ICutler et al] (l2003l) calculated the "rigidity parameter" b 
for a realistic NS structure, with a solid crust afloat on a liq- 
uid core. They solved the strain field that develops as the NS 
spins down and found that b ~ 10"^, two orders of magni- 
tude below the result found by Bavm & Pines ( 1971 1) for asim- 
plified model. On the other hand, Horo witz & Kadaul (|2009|) 
recently found through N-body simulations that the Coulomb 
lattice of the NS crust can support a maximum strain angle 
6c ~ 10', three orders of magnitude high er than the value es- 
timated bv lSmoIuchowski & Welch! (1 19701) . Thus, for an NS of 
~ 1 AMq, the time-averaged crust-cracking luminosity is L^c ~ 

1026p_2o/PLs erg s ' 
sured in units of 10" 
liseconds. For a representative classical pulsar (PSR B0950-I-08), 
P ~ 250 ms and P ~ 10"'^ s s^', so L^c ~ 10^^ erg s^', while for 
the MSP J0437-4715, P = 5.76 ms and P = 5.73x10^20 s s"', so 
Lcc ~ 4 X 10^^ erg s"' . Comparing these results with the thermal 
emission from a pulsar with Ts ~ 10^ K, L ~ 10^^ erg s"' (i.e., 
potentially detectable thermal emission from an NS in the so- 
lar neighborhood by the Hubble Space Telescope), we conclude 
that the crust-cracking mechanism does not produce detectable 
heating. 

Additionally, the high critical strain angle obtained by 
iHorowitz & Kadaul (l2009l) requires the star to have an initial 
deformation eo > 9c ~ 10 ' in order to cause any cracking 
of the crust. However, for plausible initial rotation periods of 
classical pulsars (Pq > 15 ms), their initial deformation is only 
eg ~ P\IPI) < 10"^, where Pk ~ 0.5 ms is the Keplerian period 
of the star Hence, the crust-cracking mechanism is never acti- 
vated in classical pulsars and probably operates only in MSPs 
with Pq <2 ms. If the latter were the case, and now considering 
that all stresses in the crust are suddenly released, the internal 
thermal energy of the star is increased by ~ BO^ ~ lO"*** erg, 
corresponding to an internal temperature of lO^K and a surface 
temperature of ~ 5 x 10^ K (Potekhin et al ll997h . which is dis- 
sipated within less than 10^ yr. Hence, in MSPs of 10^"^ yr, the 
increase in temperature due to catastrophic cracking is given by 
a narrow peak in the thermal evolution, which is unlikely to be 
detected because of the short timescale involved. 



2.4. Vortex creep 

The relatively low temperatures in the interior of NSs induce 
the formation of neutron Cooper pairs. These form a condensate 
with a macroscopic wave function. A consequence of this is that 
the vorticity in the superfluid must be concentrated in discrete 
vortex lines, whose microscopic distribution allows the super- 
fluid to approximate a macroscopic rigid rotation. 

As the star spins down, the vortex lines must move outward. 
As they move through the inner crust, they are pinned to the nu- 



Table 1. Vorte x-nuclei interaction pararn e 
classical model ("Ponat i & Pizzocherdl2004) 
potential (a) and Gogny potential (b). 



arameters in the semi- 
with the Argonne 



Zone p [g cm '] Rws [fm] ^ [fm] 



EZ [MeV] E^^l [MeV] 






1 


1.5xl0'2 


44.0 


6.54 


- 


2 


9.6 x 10'- 


35.5 


7.25 


- 


3 


3.4 x 10'3 


27.0 


8.54 


5.2 


4 


7.8 X 10'3 


19.4 


11.71 


5.1 


5 


1.3 X 10'* 


13.8 


8.62 


04 



7.5 
5.9 



Notes, p is the density of each zone, R^s is the radius of the Wigner- 
Seitz cell, f is the vortex coherence length, and E^p is the vortex-nuclei 
pinning energy. 



clear lattice until a critical velocity difference between the su- 
perfluid and the crust is reached. In this process, the pinning 
and unpinning of the vortex lines with respect to the nuclei of 
the crystal lattice release en ergy that heats th e star. The energy- 
dissipation rate is given by (lAlparetal.ll 19841) 



L,c = mi 



(8) 



where J ^ chip , with Ij, the moment of inertia of the pinning 
layer and lOc- - (fi.s - ^c)cr an average over the pinning zone of 
the critical lag between the angular velocity Qc of the crust and 
the superfluid rotation rate Q^. 

iDonati & Pizzocherol (|2004|) calculated the vortex-nucleus 
interaction in the inner crust of NSs with a semi-classical model. 
The density-dependent neutron pairing gaps used in the cal- 
culations are obtained from the Argonne potential and Gogny 
effective interaction. Table [T] shows the pinning energy calcu- 
lated with this model f or five zones of the inner crust. Similarly, 
lAvogadro et al.l (l2008l) calculated the vortex-nucleus interaction 
in the inner crust based on a Hartree-Fock-Bogoliubov quantum 
mean field theory. Table|2]shows the pinning energy for this ap- 
proach. An important result of this, contrary to the prediction of 
the previous model, is that pinning of vortices on nuclei is fa- 
vored at low density in the inner crust. We used both results to 
calculate the excess angular momentum J, which determines the 
luminosity for the vortex creep mechanism. In addition, we used 
Eq. (58) of lAlpar et all 019841) in the limit E^p » kT, 



_^7T r E 

3 JpK. 



EnpP_ 

mws 



dr 



(9) 



where Ej^p is the pinning energy of a vortex on a nucleus, r is the 
radial coordinate, k is the quantum of circulation of each vortex, 
^ is the vortex coherence length, and Rws is the radius of the 
Wigner-Seitz cell. In order to calculate this integral, we generate 
NS structure models for specific masses and linearly interpo- 
late the values of Table [1] and |2] according to these models. For 
simplicity, we do not take relativistic effects into account be- 
cause the corrections involved are minor, i.e. (1 - rg//?)'^^ ~ 0.8, 
with Tg the Schwarzschild radius. In this way, for an NS of M ~ 
1 .4Mq and a typical range of equations of state we find that the 
excess of angular momentum is 7 ~ (lO'^-'-lO"'^) ergs. Thus, the 
vortex-creep luminosity L^c - (10^^ - 10^')|Q_i4| erg s"', where 
Q_i4 is the angular velocity derivative in units of 10"'"* s"^. This 
is similar to the luminosity inferred from the observation of the 
MSPJ0437-4715. 
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Table 2 . Vortex-nuclei interact ion parameters in the quantum ap- 
proach ("Avog adro et alJl2008h . with the SLy4 potential (a) and 
Skm* potential (b). 

Zone p [g cm-3] R„s [fm] ^ [fm] 4"^ [MeV] 4*^ [MeV] 



10 



1.9 X 10'2 
3.3 X 10'2 
6.6 X 10'2 
1.3 X 10" 



42.9 
40.3 
37.2 
33.2 



6.63 
6.85 
7.11 
7.60 



1.08 
1.20 



1.51 


^ 


3.85 


^ 


1.63 


o 



Notes. Variables are as explained for Table [T] The v alues of ^ are inter- 
polated from Table 8 of lDonati & Pizzocherol ( 120041) . 



2.5. Rotochemical heating 

In chemical equilibrium, in an NS composed of neutrons (n), 
protons ip), and leptons (I: electrons and muons), the chemical 
potentials satisfy ri„pi = fj.„ - fj.p - fi/ = 0. However, if the rota- 
tion of the star is slowing down, the centrifugal force is reduced, 
the central density of the star increases, and the chemical poten- 
tials are imbalanced, ri„pi + 0. As the equilibrium composition is 
altered, the NS will relax to the new chemical equilibrium (via 
beta and inverse beta decays), releasing energy in the form of 
neutrinos, which leave the star, and heat, which is later radiated 
as ph otons. The evolution of the chem ical imbalances is of the 
form dFernandez & Reiseneg^erll2005h 



rinpi = -A(ri„pi, T) - R„piQ.Q, 



(10) 



where the function A quantifies the effect of reactions toward 
restoring chemical equilibrium, and the scalar R„pi quantifies the 
departure from equilibrium due to the change in the angular ve- 
locity Q of the star 

The luminosity generated by this mechanism is L = r?;,,^/, 
where T = r„^p/v - Tpi^ny. Here, r„^piy is the rate of reac- 
tions (integrated over the core) that transform the neutrons to 
protons and leptons through direct or modified Urea reactions, 
and Tpi^„y is the rate for the opposite process. In this way, the 
evolution of the internal temperature with rotochemical heating 
is given by the solution of the coupled differential Eqs. ([1]) and 



iReiseneggeiJ (Il995 f) found that if the angular veloc- 
ity Q varies slowly over the time required to cool the 
star and achieve chemical equilibrium, the star reaches a 
quasi-steady state, where heating and cooling are balanced. 
[Fernandez & Reiseneggeil |2005) calculated the simultaneous 
solution of t - fjupi - for a typical range of equations of 
state and found that in an NS with a non-superfluid core and 
with modified Urea reactions, the photon luminosity in the quasi- 
steady state depends only on the period and its derivative. 



L'^ ^ (10''' - 10'') 



P 



20 



8/7 



ergs 



(11) 



which is close to matching the observation of MSP J0437-4715. 
Here, P_20 is the period derivative measured in units of 10"^" and 
f tns is the period in milliseconds. The characteristic timescale to 
reach this quasi -steady state is Te,„ ^ 2x 10^(f ^ /P_7,n)^^^ yr On 
the other hand, iPetrovich & Reiseneggeil (1201 Oh considered the 
effects of nucleon superfluidity on rotochemical heating. They 
found that the chemical imbalances grow up to a value close 
to the energy gaps, which is higher than in the nonsuperfluid 
case. Therefore, the surface temperatures predicted with Cooper 



- 


1 1 1 1 1 
Vortex creep 


■ 


mod. Urea 




: 


. dir. Urca~~-- 


\^ 


■ 


"\ "^'^^^'^^c:;^ 


^^^^^"^="^^^^^^^^5^ 






^^^^ 




1 >■■''■:. . 1 


■ 



log t [yr] 
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Fig. 1. Thermal evolution with the vortex creep mechanism. All 
curves show the interior temperature as a function of time, for 
stars with mass M = I.4M0, initial temperature T - 10" K, and 
initial period Pq - 1 ms. The magnetic fields B = 10^ G and 
B = 10" G correspond to MSPs and classical pulsars, respec- 
tively. The excess angular momentum used is J - 10'*-' erg s. 
The solid lines show the evolution using the AI8-1- 6v + UIX* 
EOS, with only modified Urea reactions, and the dashed lines 
show the evolution using BPAL 22 EOS, with direct Urea reac- 
tions. The dotted lines show the evolution with passive cooling 
for direct and modified Urea reactions. 



pairing are higher. For simplicity, we here only consider the non- 
superfluid case. However, we must not lose sight of the latter 
result. 



3. Effects 

In order to solve the thermal balance equation, Eq. ([T), 
and generate the evolutionary curves, we use the code of 
[Fernandez & Reisene^gej ( l2005h . This considers realistic equa- 
tions state (EOSs), with a conventional NS divided into two re- 
gions: an isothermal interior and a thin envelope. In the interior, 
it considers a core composed of neutrons, protons, electrons, and 
muons, but ignores the potential Cooper pairing effect. The neu- 
trino emissivity is generated by modified Urea reactions and di- 
rect Urea reactions, neglecting any other neutrino emission pro- 
cesses. In order to model the envelope, it uses the relation be- 
tween internal and surface temperature f rom the fully accreted 
envelope model of iPotekhin et alJ (Il997h . The rotational evolu- 
tion is assumed to be due to magnetic dipole radiation, without 
magnetic field decay. 

3.1. Vortex creep 

Figure [T] shows the resulting thermal evolution for the vortex 
creep mechanism for representative parameters of both classi- 
cal pulsars and MSPs, with one EOS allo wing only for mod- 
ified Urea reactions (A18-H(yv-HUIX* EOS. lAkmal et al.lll998h 
and another one a llowing for direct Urea reactions (BPAL 22, 
iPrakash et alJll988i) . As the residual spin-down at old ages is 
greater in pulsars with weaker magnetic fields, the predicted 
temperature is higher in the MSPs than in the classical pulsars. 
At very late times (f > 10^ yr), the cooling is dominated by 
photon emission, which balances the heat generation by vortex 
friction. Therefore, the evolution of temperature becomes inde- 
pendent of the previous thermal history and the type of Urea 
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Fig. 2. Thermal evolution with rotochemical heating. All curves 
correspond to stars with mass M - IAMq, A18 + 6v + UIX* 
EOS (which allows only for modified Urea reactions), initial 
temperature T = 10^ K, and magnetic field B = 10^ G. The 
solid lines show the evolution of the internal temperature T and 
the chemical imbalance 77,,;,^ , and the dotted line shows the pas- 
sive cooling. 



10 



Rotochemical heating (Classical pulsars) 

__77,^(P„=lms) 




Fig. 3. Thermal evolution with rotochemical heating. All curves 
coiTespond to stars with mass M - IAMq, A18 + Sv + UIX* 
EOS, initial temperature T = 10" K, and magnetic field B - 
10" G. The solid and long-dashed lines show the evolution of 
the temperature and the chemical imbalance with initial periods 
Po = 1 ms and Pq - 15 ms, respectively. The dotted line shows 
the passive cooling. 



reaction, and only depends on the current value of Q, with the 
photon luminosity L = AnR^crT^ - J\Cl\. The rotational evolu- 
tion by magnetic dipole radiation yields Q. oc f ''^ in the limit 
Q <s: Qi). Thus, the asymptotic decrease of the surface tempera- 
ture with time in classical pulsars and MSPs is T, oc r^^^. 



3.2. Rotochemical heating 

The thermal evolution for MSPs with rotochemical heating and 
modified Urea reactions is shown in Fig.|2] For this mechanism, 
the thermal evolution is coupled with the evolution of the chem- 
ical imbalances. Because the magnetic fields in MSPs are rela- 
tively weak, B ~ 10** G, the chemical imbalances induced by 
the spin-down grow slowly, causing chemical reactions at high 



ages, f ^ 10^ yr An important prediction of this mechanism is 
that in the MSP regime, the star arrives at a quasi-steady state, 
where heating and cooling balanced, so the thermal evolution is 
independent of initial conditi ons and only depends on the cur - 
rent value of the product D.Cl ([Fernandez & Reiseneggeij[2005h . 
When direct Urea reactions are present, the evolution is qual- 
itatively similar, but the temperatures are strongly reduced, as 
shown in Sect. |4] 

The evolution of clas sical pulsars is shown in Fig.|3](see also 
iReisenegger et al.ll2007l) . As in the previous figure, only modi- 
fied Urea processes are active. Owing to the relatively strong 
magnetic fields present in this regime, B ^ 10" G, the ro- 
tational energy of the star is quickly consumed in the begin- 
ning of the thermal evolution by the magnetic dipole radiation. 
Because of this, high chemical imbalances are built up, propor- 
tional to the amount of rotational energy lost. Thus, for faster 
initial rotation, the chemical imbalance at later times will be 
higher. At ages above ~ 10'^ yr, the right-hand side of Eq. (fTOl l is 
dominated by the reactions restoring the chemical equilibrium. 
When only modified Urea reactions are present and 77 » kT, the 
chemical imbalance evo lves according to ?) oc -F oc -rf (see 
iFernandez & Reisene^^ej l2005h . yielding 77 oc r '/^. Because 
in the last epoch of the thermal evolution the cooling of the 
star is caused mainly by photon emission, the luminosity is 
L oc T^ oc Trj oc rf . Thus, the surface temperature decreases 
with age as Ts oc r^^^, only slightly more slowly than with vor- 
tex creep. 

4. Comparison with observations 

In order to identify the thermal emission at T ~ 10^ K from the 
whole NS surface, it is necessary to obtain ob servations in the 
optical and specially in the ultraviolet range (Kar galtsev et al.l 
2004; Zavlin & Pavlov 2004). Thus, to confront the theoretical 
cooling models of old NSs with observations, we select the 
best and most restrictive observations in these bands, of seven 
very old pulsars: three MSPs and four classical pulsars. Among 
the MSP s, only J0437-4715 has detected thermal emission 
(Karg alts evetal.1 [2004), and J2124-3358 ( Mign ani & Beckgi 
12004 and J0030-I-0451 dKoptsevich et all l2003h have good 
upper limits on their temp eratures. The four selected clas- 
sical pulsars, B 1 929-H 10 dBecker et alj 12006)) . B0950+08 
(Zavlin & Pavlov '2004; Zharikov et al. 2004), Bl 133+16 
(Zharikov etal. 2008) . and J0108-1431 (Mig nani eTaP l2003h 
only have upper limits on the thermal emission from the whole 
NS surface. 

The upper panel of Fig. |4] shows the thermal evolution of 
MSPs including vortex creep and rotochemical heating with 
modified Urea reactions. For the vortex creep mechanism the pa- 
rameter J is highly uncertain, so we adjust it to the lowest value 
that is consistent with the thermal emission from MSP J0437- 
4715. This constrains it to 7 > 5.5 x 10'*"' erg s if vortex creep 
is postulated as the main heating mechanism for this pulsar For 
rotochemical heating, the surface temperature of this MSP can 
be quite precisely p redicted, considering its mass M - \ .76Mq 
(IVerb iest et al' 2008) and an interior model given by the A18 + 
6uj¥ UIX * EOS. This pre diction is 1 .Icr below the observation 
of Karg altsev et all (l2004 . 

Similarly, the bottom panel shows the thermal evolution of 
the classical pulsars. For the vortex creep mechanism, we as- 
sumed that the parameter J is universal. Thus, to solve the ther- 
mal evolution in this regime, we used the previous value of J im- 
posed by the MSP J0437-4715, resulting in a "phenomenologi- 
cal" temperature prediction slightly above the constraints for the 
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Fig. 4. Evolution of the surface temperature for a neutron 
star with vortex creep (long-dashed lines), rotochemical heat- 
ing (solid lines), and passive cooling (dotted lines). All curves 
correspond to stars with the K\% + 6v + UIX* EOS and mod- 
ified Urea reactions. The error bar shows the temperature mea- 
sured for the MSP J0437-4715 and the arrows show the upper 
limits on the thermal emission for specific pulsars. Top panel: 
The curves correspond to MSPs with mass M - \ .16Mq, mag- 
netic field B = 2.8 X 10*^ G, and initial temperature T = lO*^ K. 
Bottom panel: The curves correspond to classical pulsars with 
mass M - \AMq, magnetic field B = 2.5 x 10" G, and initial 
temperature T - 10" K. The abscissa corresponds to the spin- 
down time (fj - Q/2|Q|). The initial periods for rotochemical 
heating are labeled on each curve. 
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Fig. 5. Evolution of the surface temperature for a neutron 
star with vortex creep (long-dashed lines), rotochemical heat- 
ing (solid lines), and passive cooling (dotted lines). All curves 
correspond to stars with the BPAL 32 EOS and direct Urea re- 
actions. The error bar shows the temperature measured for the 
MSP J0437-4715 and the arrows show the upper limits on the 
thermal emission for specific pulsars. Top panel: The curves 
correspond to MSPs with mass M - 1.76Mq, magnetic field 
B = 2.8 X 10*^ G, and initial temperature T = 10** K. Bottom 
panel: The evolutionary curves correspond to stars with mass 
M - IAMq, magnetic field B = 2.5 x 10" G, and initial tem- 
perature T - 10" K. The abscissa corresponds to the spin-down 
time {ts = Q/2|Q|). The initial period for rotochemical heating 
is P() = 1 ms. 



pulsars B1929-I-10 and B0950-I-08. On the other hand, as we saw 
in Sect. 13.2! the evolution with rotochemical heating depends 
on the initial period of rotation, Pq, for NSs with relatively high 
magnetic fields. Because of this, rotochemical heating can easily 
accommodate substantially lower temperatures if more restric- 
tive observations of classical pulsars are available. The evolution 
at late times (t > 10** yr) is remarkably similar for both mecha- 
nisms. We showed that the relation between surface temperature 
and age is Ts cc r^/^ for vortex creep (Sect. 13.1b . very slightly 
steeper than the relation Ts oc f '^^ for rotochemical heating with 
modified Urea reactions (Sect. lTTb in the latter stage of the ther- 
mal evolution. 

Figure|5]shows the thermal evolution considering the fast di- 
rect Urea processes. For the vortex creep mechanism, we used 
the same excess of angular momentum J as in the previous evo- 
lutionary curves, fixed to be consistent with the observed temper- 
ature of MSP J0437-4715. Because for ages S lO'^yr this mech- 
anism does not depend on the Urea process type, the tempera- 
tures of the MSPs (upper panel) are the same as with modified 



Urea reactions. However, in classical pulsars (bottom panel) the 
thermal evolution is sensitive to the direct Urea processes and 
generates a predicted temperature slightly lower than the limit 
for PSR B0950-H08. In this way, a more sensitive observation of 
this pulsar could discard vortex creep as the main source of the 
thermal emission detected from the MSP J0437-4715. On the 
contrary, the temperatures generated by rotochemical heating, 
both for classical pulsars and MSPs, are strongly reduced if di- 
rect Urea processes are considered. The temperature measured in 
J0437-4715 requires that the neutrino emission is produced only 
by modified Ure a react ions (unless substantial Cooper pairing 
gaps are present: fPetrov ich & Reisenegger 2010) if rotochemi- 
cal heating is the main source of heat. 

Finally, we compare the excess angular momentum, J, that 
we estimated from the pinning energies of Tables [1] and |2] 
with the observations. In order to do this, we numerically in- 
tegrated Eq.|9]over the inner crust of an NS with canonical mass 
M = 1.4M0 and the A18-H(?v-h UIX* EOS. In this way, wit h 
the pinning energy 



He Ais-HOv-H UiA'^ fcuti. in this way, wit h 
estimated by iDonati & Pizzocherd (l2004l) . 
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we obtain J - 3.8 x lO**^ erg s for the Argonne interaction and 
J - 5.9 X 10^"* erg s for the Gogny interaction. Similarly, but with 
the pinning energy estimated by Avogadro et al. (2008), / - 
1.2 X 10'*^ erg s with the SLy4 interaction, and 7 = 6.8 x 10"*^ erg 
s with the Skm* interaction. Comparing these results with the 
value J - 5.5 X 10"*^ erg s, which is marginally compatible 
with both the observed temperatures of MSP J0437-4715 and 
the upper limits for the old classical pulsars, we fin d that the 
semi-classical model of iDonati & Pizzocherol ( |2004|) overesti- 
mates the value of t he excess angular mom etum J, while the 
quantum approach of lAvogadro et al.l (l2008h . with the SLy4 in- 
teraction, roughly agrees with the value inferred from observa- 
tions. 



5. Conclusions 

We studied five heating mechanisms that can be operating in old 
neutron stars: magnetic field decay, dark matter accretion, crust 
cracking, vortex creep, and rotochemical heating, and compared 
them with pulsar observations. 

We found that magnetic field decay, dark matter accretion, 
and crust cracking cannot produce detectabl e heating in old pul- 
sars. Owing to the high yield strain angle (iHorowitz & Kadaul 
[2009), the crust cracking mechanism does not operate in clas- 
sical pulsars, and probably only operates in MSPs. The vortex 
creep and rotochemical heating can be important both for classi- 
cal pulsars and MSPs. 

In the evolutionary curves with vortex creep, and with the 
excess angular momentum, J, set to the lowest temperature com- 
patible with the thermal emission detected from MSP J0437- 
4715, the predicted temperature turns out to be very near the ob- 
servational upper limits for the classical pulsars B1929+10 and 
Bl 133-1-16. Likewise, rotochemical heating with modified Urea 
reactions and no Cooper pairing is only 1 .7 cr below the tem- 
perature measured in the MSP J0437-4715. However, the tem- 
perature prediction of rotochemical heating can be raised if a 
superfluid core is considered in the model. 

The prediction of pinning energies in the inner crust of NSs 
via the semi-classical mode l for the vortex-nuclei interaction 
dPonati & Pizzocherol |2004 |) overestimates the temperatures of 
B0950-I-08 and B 1929H-10. The recent estimations of pinning en- 
ergies via a quantum approach (Avogadro et al. 2008) are consis- 
tent with the observations of classical pulsars and MSPs. Finally, 
more stringent constraints on the temperature of some classical 
pulsars such as B0950+08 could rule out the vortex creep mech- 
anism as the main source of the thermal emission detected in the 
MSPJ0437-4715. 
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